Volt age- driven quantum oscillations of conductance in graphene 
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Abstract 

Locally-gated single-layer graphene sheets have unusual discrete energy states inside the poten- 
cy : 

tial barrier induced by a finite-width gate. These states are localized outside the Dirac cone of 

continuum states and are responsible for novel quantum transport phenomena. Specifically, the 
OO ' 

' longitudinal (along the barrier) conductance exhibits oscillations as a function of barrier height 

O ' 

and/or width, which are both controlled by a nearby gate. The origin of these oscillations can 



be traced back to singularities in the density of localized states. These graphene conductance- 



m 
o 

oscillations resemble the Shubnikov-de-Haas (SdH) magneto-oscillations; however, here these are 



driven by an electric field instead of a magnetic field. 

X" 

03 1 PACS numbers: 73.63.-b 
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I. INTRODUCTION 



The unusual and rather remarkable transport properties of graphene continue to at- 
tract considerable attention. Soon after its experimental discovery^, studies found: un- 
conventional quantum Hall effect^; the possibility of testing the Klein paradox-; specular 
Andreev reflection and Josephson effect^; new electric field effects^; intriguing electron 
lensing^; and other fascinating phenomena (see, e.g., recent papers^^^^ 1 ^^^ and 
references therein). Studies of graphene are also inspired by their potential application in 
nano-electronic devices, since an applied electric field can vary considerably the electron 
concentration and have both electrons and holes as charge carriers with high mobility. 

The subject of the present study, which is a logical continuation of recent work^, is an un- 
usual novel transport effect, namely, voltage- driven quantum oscillations in the conductance 
of a single-layer gated graphene. These oscillations originate from a new type of electron 
states in graphene. When a graphene sheet is subject to nearby gates, these create an energy 
barrier for propagating electrons. Here we explicitly demonstrate that, in contrast to non- 
relativistic quantum mechanics, where localized states can exist only inside quantum wells, 
Dirac-like relativistic electrons in graphene allow energy states localized within the barrier. 
We show that the energy e(q y ) of the localized states (versus the wave vector component q y 
along the barrier) becomes non-monotonic if 

VqD > irhvp, 

where Vq and D are the barrier height and width correspondingly, and vf is the Fermi 
velocity. This produces singularities in the density of localized states for energies where 

de/dq y = 0. 

When the magnitude and/or width of the barrier changes, the locations of the singular- 
ities move and periodically cross the Fermi level, generating quantum oscillations in the 
longitudinal (along the barrier) conductance as well as in the thermodynamic properties of 
graphene. This situation resembles the well known physical mechanism for Shubnikov-de- 
Haas (SdH) magneto-oscillations (see, e.g., Ref.— »^). Indeed, electrons in the conduction 
band of a 3D metal subject to a strong magnetic field have an equidistant discrete energy 
levels (Landau levels) separated by the cyclotron energy. The corresponding density of states 
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has singularities at the Landau levels. When the magnetic field is changed, the positions of 
the Landau levels move and pass periodically through the Fermi energy. As a result of this, 
the population of electrons at the Fermi level also changes periodically, giving rise to the 
quantum oscillations of both the transport and thermodynamic properties of a metal. One 
should notice, however, a few important differences. First, in the context of gated graphene, 
oscillations are induced by the electric field, while the corresponding SdH oscillations are 
driven by a magnetic field. Second, localized energy states in graphene are non-equidistant 
and the resulting density of states has a rather complicated energy dependence. Thus, the 
corresponding oscillations in the conductance inherit all these unusual peculiarities. 



II. LOCALIZED ENERGY STATES IN A BARRIER 

The tunneling of relativistic particles in graphene across a finite-width potential barrier, 
and its corresponding conductance, has been recently studied (see, e.g., Refs.-^ 1 ^ 1 ^ 1 ^) . Here 
we consider another conduction problem, namely, electron waves that propagate strictly 
along the barrier and damp away from it. More specifically, we consider electron states in 
graphene with a potential barrier located in a single-layer graphene occupying the xy-plane 
(see Fig. [TJ. For simplicity, we assume that the barrier V(x) has sharp edges, 
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V(x) = (1) 



Electrons in monolayer graphene obey the Dirac-like equation (hereafter H = 1), 

= Hiff, H = -iv F a ■ V + V{x) , (2) 

where <x = (a x , a y ) are the Pauli matrices. We then seek stationary spinor solutions of the 
form, 

i/j(x, y) = ifj(x) exp(-iet + iq y y) , (3) 

with energy e and momentum q y along the barrier. We focus on states with \q y \ > |ac| = 
\e\/vp- In this case, the electron waves satisfying Eq. (T5]) damp away from the barrier, and 
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FIG. 1: (Color online) Schematic top view of a graphene sheet (yellow rectangle) placed under 
voltage gates indicated by the grey block rectangles. Bottom: gate-induced potential energy barrier 
V(x) in graphene. 

the components ip\ and ip2 of the Dirac spinor can be written in the from 

a exp [k x (x + D/2j\, x<-D/2, 
b exp (iq x x) 

+ c exp (—iq x x) , \x\ < D/2 , 
k d exp[- k x (x- D/2)] , x>D/2, 



^i(x) 



(4) 



' a ^f- exp [k x (x + D/2)] , x < -D/2 



(5) 



c+Qy 

—b exp(iq x x + id) 

+ c exp(— — i9) , |x| < .D/2 

-d^exp[-k x (x-D/2)], x>D/2, 

with real /c x = {q 2 — k 2 ) ^ 2 and g x = ((« — V/D) 2 — q 2 ) 1 ^ 2 ■ Here V = VqD/v f is the effective 
barrier strength and tan 6* = q y /q x . 

Matching the functions ipi(x) and ^(x) at the points x = ±D/2, we obtain a set of 
four linear homogeneous algebraic equations for the constants a, b, c, and d. Equating the 
determinant of this set to zero, we obtain a dispersion relation for the localized electron 
energy states, 



F(e, q y ) = tsna(q x D) 



k x q x 



«( V /D-«) + g2 



0. 



(6) 



The spectrum of localized states in graphene [Eq. is shown by the solid black curves 
in Fig. [21 for dimensionless variables 

Q = q y D, £ = eD/v f . (7) 

This spectrum consists of an infinite number of branches S n (Q). Each of these branches 
starts from the lines £ = ±|<5| ( re d solid straight lines in Fig. [2]) at 

8 = V/2 - vrV/2V (8) 

and tends asymptotically to the line 

£ = V-Q 

with increasing Q (dashed red line in Fig. [2]). Furthermore, a particular branch of the 
spectrum starts at the point (Q = 0, £ = 0) and also tends to the line £ = V — Q. 

The behavior of different branches of the spectrum depends on the barrier strength V. 
If V < 7r/2, then all branches satisfy £ < 0. Localized states with positive energies appear 
only for V > ir/2. When V increases, new branches in the spectrum with positive energies 
appear. When V is within the interval 

(n + 1/2)tt < V < (n + 3/2)tt, 

the number of branches with £ > is n + 1, for n = 1, 2, 3, . . .. It is worth emphasizing that 
each of the branches with positive energy has a maximum £™ ax at a certain wave number 
Q = Qn ax - Near these points the group velocity of localized electron waves tends to zero, 
which resembles the stop-light phenomena found in various medial. The localized states 
can also be observed in graphene when a voltage is applied to produce a potential well^. 

Note that defect-induced localized electron states in graphene and the enhancement of 
conductivity due to an increase of the electron density of states localized near the graphene 
edges were recently reported^. In contrast to these examples, the electron states studied 
here are localized within the barrier and also these are tunable, i.e., the energy levels can be 
shifted by charging the barrier strength (e.g., via tuning a gate voltage). 
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III. DENSITY OF LOCALIZED STATES 



To calculate the density of electron states p(e), we use the general formula p{e) = J2 a ^( £ ~ 
e a ), where the index a labels the quantum state and S(x) is Dirac's delta- function. Using 



£(. . .) = 4L x L y (27r)- 2 J dk x dk y (. . .) 



(9) 



for a continuum spectrum one finds the already familiar expression 



Pcont(£) = P |£|> Pa 



(10) 



TTVpD ' 

where L x and L y are the lengths of the graphene sheet in the x and y directions, respectively. 
For localized energy states, we obtain: 



D 



P 



JS)=2p -J2 



dS n (Q) 



dQ 



(11) 



Sn(Q)=£ 




FIG. 2: (Color online) Electron energy spectrum in graphene obtained for positive Q and V = 16. 
The Dirac sea of delocalized states (continuum spectrum) is marked by the grey regions. The 
branches of the spectrum for localized states are shown by solid black curves between the straight 
solid and dashed red lines. There are no states in the forbidden (white) regions. The increase or 
decrease (schematically indicated by the upward and downward vertical blue arrows) of the Fermi 
level £p (marked by the horizontal dashed blue line) results in a periodic change in both the density 
of states at £ = £p, and also in the conductance. 
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FIG. 3: Dimensionless oscillating parts of the density of states p loc /p at the Fermi level (upper 
panel) and conductance g loc /g cont (lower panel) versus the strength V of the potential barrier for 
D/L x = 0.1, £p = 1 (main panels) and £ f = 5 (insets). The total conductance is g = g cont + g loc . 



where n runs over the positive roots of the equation £ (Q) = £. The function p loc (£) exhibits 
two types of peculiarities. First, increasing £, the jumps or steps (each of magnitude 2D / L x ) 
in p loc {£)/p occur at the points, given by Eq. ([8]), where new branches of the spectrum arise 
or disappear. More importantly, singularities are observed when £ = £™ ax , where |d£/d<3| _1 
in Eq. ffTT]) diverges. 

The locations of the singularities shift when changing the barrier strength V. Therefore, 
they periodically cross the Fermi level £p. This produces quantum oscillations in the density 
of states at the Fermi energy, which are seen in the upper panel of Fig. [31 showing p loc (£)/ p 
versus the effective barrier strength V. 
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IV. KUBO FORMULA AND CONDUCTANCE 

A. Kubo expression for the conductance in graphene 

When studying transport, within linear response theory, one usually starts from the 
current-response function, 

J^(x,x') = -z^-OTr{^jJ(x),jf(x')]} , (12) 

where x = (r, t), ~&(t) is the Heaviside step-function, g is the equilibrium density matrix, 
and 

jjf (r, t) = exp (i#i)j M (r) exp (-iHt) 

is the current operator in the Heisenberg representation with the Hamiltonian taken from 
Eq. (j2J), and where [...,...] stands for the commutator. For electrons with a linear Dirac 
spectrum, one finds 

j^r) = ev F ft(r)a^(r). (13) 
Equation f|T2|) is used to define the frequency-dependent linear conductance as 

= ^^TT // dr dr '^( r > r » • (14) 

Here 3ft stands for the real part of a complex number. 

We expand the fermionic field operator i/j(y, t) in terms of exact eigenfunctions [Eq. (j3J)], 
namely, 

^(r, t) = ^ V«( r ) exp (-ie a ta a ), (15) 



and then perform quantum averaging in Eq. ( 1121) with the help of Wick's theorem and the 
relation Trj^a^ag} = 5 aj g/(e a ), where 

/(e) = l/[exp[(e-e F )/r] + l] 

is the Fermi occupation function. Performing a Fourier transform and using 

Vl[i/(e-e' + u + iQ)] = ir5(e-e f + w) (16) 

Eq. (|T4l) . reduces to 

vr(e^) 2 /■+« /(e+)-/(e_) 



^(w) = / de 



xTr [a^{e + - H) rr ,aJ(e„ - H) r , r } , (17) 
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where e± = s ± uo/2 and the trace incorporates spatial integrations. The operator delta- 
functions can be directly related to the single-particle Green's functions 

G e (r,r') = (r\(e-H)- 1 \r') 

according to 

S(e - H)rr> = ^ [G a e(r, r') - G r £ (r, r')] , (18) 

where the superscript a/r stands for the advanced/retarded component, respectively. As a 
result, one finds for the conductance 

^.Mr^ (19) 



x 



Tr{a, [<%(r, r') - GI + (r, r')] a v [G r £ _(r', r) -G a £ _(r', r)]} . 



Next we incorporate disorder by introducing the one-particle scattering time r, for Dirac 
fermions, into the Green's function, 

(G^ a ) dis ^(6-H± i /T)'\ (20) 

which enters through the imaginary-part of the corresponding self-energy. The subindex 
"dis" refers to disorder. Furthermore, we factorize the average of the product of two Green's 
functions by the product of their averages, 

(G; + G a £ _) dis ^(G r £+ ) dis (G a £ _) dis . (21) 

This assumption should be valid for weak disorder and together with Eq. (|2"U|) is equivalent 
to the self-consistent Born approximation. 



B. Conductance along the barrier 

We now focus on the along-the-barrier (fi = v = y) conductance for the geometry shown 
in Fig. [TJ At zero temperature, T — > 0, when f(e) = $(ef — s) and the e integration is 
bounded by the frequency o>, for the average dc-conductance g = (g yy (uj — > 0))di S we find 
(per spin and per valley): 

9 = 9 con t + 9 loc ■ (22) 
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The first contribution g cont here comes from the extended electron energy states with corre- 
sponding density of states taken from Eq. ffTU]) . and reads explicitly (now keeping h) as 

7re 2 L 



9 C 



1QHL 



y l 



e f t H — 1 — e f t arctan 

71 V E F T 



(23) 



At the neutrality point, Ef = 0, from Eq. fl23l) one recovers a universal (i.e., scattering time 
r-independent) result g cont = a miri (L x /Ly), where a m \ n = (n/8)(e 2 /h) is the minimal conduc- 
tivity, which received considerable attention in a number of recent studies (e.g., Refs.— 
Away from the neutrality point, the conductance growths linearly with the Fermi energy, 
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{ire 2 /16h)(L x /L y )e F r. 



(24) 



The novel result of the present study is the oscillatory part g loc , which originates from the 
electron states localized within the barrier. It can be expressed, with the help of Eq. (fTTI) . 
as follows: 



9 loc 



2e 2 D 



POO 

/ d£ 


dQ 




Jo 







M(S)t] 2 



2 ' 



(25) 



where 



M{£) 



dx 

u 



is the matrix element constructed from the wave- functions of localized states, Eq. (jl])-©, 
and f] = D jvpT. The remaining integration in Eq. (|25]) is simplified realizing that everywhere 
away from the integrable square-root singularities of |dQ/d£ n |, the ^-dependent function is 
peaked at the Fermi energy, whereas M(S) is smooth. Thus, one finally finds, 

dQ 

- JVI 

fi'cont L x 



d£ n 



(26) 



where the conductance g loc is normalized to its continuous part taken away from the neu- 
trality point, namely, where g cont oc te f . Note that 

D 



S f 



E F - 



V F 



The derivative entering Eq. ([26]) can be calculated with the help of the dispersion equation 
(ED as (dQ/dS) = -(dF/dS)/(dF/dQ), and reads 

V-2£ + {V-£) v / Q 2 -£ 2 



dQ 
d£ 



Q 



iy-E)E-Q 2 -Q 2 ^Q 2 -^ 



(27) 
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FIG. 4: Dimensionless oscillating parts of the density of states p loc /p at the Fermi level (inset) 
and conductance g loc /g cont (main panel) versus the Fermi energy £p, for D/L x = 0.1, and V = 16. 

The oscillatory nature of <7 loc (V, Bp) is illustrated in the lower panel of Fig. [3J The essential 
observation, which follows from Eq. (1261) . is that the longitudinal conductance traces the 
peculiarities in the density of localized states and opens a direct way for their experimental 
observation. It is also worth mentioning that close to the singularity of dQ/d£, meaning 
\S n — £p\ % rj, the conductance correction is regularized by the finite width of the ?7-Lorenzian 
under the integral of Eq. 

Varying the concentration of free particles with constant barrier strength, one can again 
observe oscillations in the density of states (see the inset of Fig. HJ). Thus, the part of the 
conductance originated from the localized states, also oscillates with the change of the Fermi 
energy (see main panel of Fig. HJ). 



V. CONCLUSIONS 



In summary, we predict a novel type of conductance oscillations in locally-gated single- 
layer graphene, which are related to the unusual electron states localized within a potential 
barrier. When the barrier height and/or width is varied, localized levels periodically cross 
the Fermi energy, inducing modulations in the density of states. The latter translates into 
unusual quantum oscillations of the conductance. These electric-field-driven quantum oscil- 
lations are similar to the Shubnikov-de-Haas oscillations which are produced in metals and 
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semiconductors when changing the external magnetic field. 
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